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Abstract 

We study the normal state of a 3-d homogeneous dipolar Fermi gas beyond the Hartree-Fock 
approximation. The correlation energy is found of the same order as the Fock energy, unusually 
strong for a Fermi-liquid system. As a result, the critical density of mechanical cohapse is smaller 
than that estimated in the Hartree-Fock approximation. With the correlation energy included, a 
new energy functional is proposed for the trapped system, and its property is explored. 
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I. INTRODUCTION 



The successful creation of KRb polar molecules has provided a new many-body system to 
explore |l| . A dipolar Fermi gas is quite different from a typical Fermi gas or an electron gas 
due to the long-range and anisotropic interaction between the dipoles. Theoretical studies 

n 

on various properties of the dipolar Fermi gas have been carried out p]. The Fermi surface 
is deformed into an elliptical shape jsl in the normal state of a homogeneous system. A 
p-wave superfiuid state was proposed 4|- So far most theoretical studies were done in the 
Hartree-Fock (HF) approximation which underestimates correlation effects. In this paper, 
we go beyond the HF approximation to obtain the correlation energy of a homogenous 
dipolar Fermi gas in the normal state, and use it in the density-functional description of 
inhomogeneous cases. 

This paper is organized as follows. In section II, the correlation energy of a homogeneous 
dipolar Fermi gas is obtained in the perturbation theory. We find that the ground state 
energy is lowered by the correlation effect, leading to a significantly smaller critical density 
of mechanical collapse than the HF result. In section HI, the expression of the correlation 
energy is used to construct the energy functional of a trapped dipolar Fermi gas in the 
local density approximation. This energy functional is shown to satisfy a virial theorem. 
The critical total number of polar molecules for mechanical collapse in the trap is obtained. 
Conclusions are given in section IV. 



II. A HOMOGENEOUS DIPOLAR FERMI GAS 



We consider a homogeneous dipolar Fermi gas with mass m and permanent electric dipole 
moment d polarized along the z-axis. The uniform system is described by the Hamiltonian 

H = Y. efc4«k + ^ ^q4i+q42-q«k2aki, (l) 

k QT^Okikz 

where ak and are the fermion annihilation and creation operators, = h^k'^ /{2m) is the 
fermion kinetic energy, = 4:7icP{3cos^ (pq — l)/3 is the dipole-dipole interaction for the 
transferred wavevector q, 0q is the angle between the 2;-axis and q, and V is the volume. 

In the HF approximation, the single-particle excitation energy and the Fermi surface 
are no longer isotropic due to the anisotropic dipole-dipole interaction j^, and the HF 
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Hamiltonian is given by 

1 

k ki,k2 

where the single particle energy is given by 



HhF = Xl^k4ak - ^ J2 ^ki-kankiTlkj, (2) 



£^k = Cfc + ^ Vk_k'nk' , 

^ k' 

n-k = 0{fi — ^k) is the occupation number, and the chemical potential fi is determined from 
the total number constraint = X^k ^k- In the HF approximation, the occupation number 
nk can be solved self consistently. To the first order in the interaction strength, the deformed 
Fermi surface is given by [sl 

m = kF + ^^^^{3cos'^-l), (3) 

where kp is the Fermi wavevector of the noninteracting Fermi gas. To the first order, this 
deformation does not change the volume inside the Fermi surface, nor the total energy. The 
HF approximation modifies the ground state energy starting only in the second order js], 

EffF h'^kp 1 1 mcfkp 

or in terms of the Fermion density n, 

EhF r 3 5 4 2\i/"^^^N2 ^1 /A\ 

= — — (67r )3n3 (677 13 ( — —Vnal^ (4) 

where n is the Fermion density in real space. 

Although the HF approximation is accurate for the ground state energy in the first order, 
in general it is insufficient for the second or higher-order ground state energy. Kohn and Lut- 
tinger js] showed that in Fermi systems with anisotropic Fermi surfaces due to interaction, 
the second-order ground state energy consists of two parts, one from the HF approximation 
and the other from the Brueckner-Goldstone (BG) formula. The BG formalism is the stan- 
dard perturbation theory which was used by Lee and Yang to obtain the second-order ground 
state energy of a dilute Fermi gas [7|. In the BG formalism, the second-order contribution 
to the ground-state energy of a dipolar Fermi gas is given by 

p(2) _ (0|//i|m)(m|gi|0) 

^BG - 2^ ^ ^ ' {^) 
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where | 0) and | m) are the ground and excited states of the noninteracting Fermi gas with 
eigenenergies Eq and Em, and Hi is the dipole-dipole interaction given by the second r.-h.-s. 
term in Eq. ([T]). This BG energy can be further written as 

X/ki/k2(l - /ki+q)(l - /k2-q), (6) 

where /k = 0{kF — k) is the fermion occupation number in the noninteracting ground state 
I 0). The k-space summations in Eq. ([6]) can be rescaled into integrals over dimensionless 
variables, 



eP(j md^kl /■ ,3 ,3 ,3 (cosVy - |)(cos2 0. 



(rxi(rx2d^y 



gyv^i^vjo V^y COS v^x2— xi— y 



V 32rfj y a;2^2;i- |xi+y|2- |x2-y|2 

xe{l - xi)e{l - a;2)^(|xi + y| - l)^^(|x2 - y| - 1), (7) 

where Xj = kj/Zc^? and y = (\/kp. Eq. ([7]) indicates that the BG energy density is propor- 
tional to d^n'/^. We obtain the coefficient by Monte Carlo integration [Sj, 

_p = _o.66^n3. (8) 

By definition, the correlation energy is the difference between the true ground state energy 
and the ground state energy in the HP approximation. Therefore in the second order, the 
correlation energy of a dipolar Fermi gas is given by Eq. (|8]). Its negative sign indicates 
that the ground state energy is lowered by the correlation effect. 

The total ground-state energy of a homogeneous dipolar Fermi gas up to the second order 
can be obtained by adding Eq. ([8]) to Eq. (jlj), 

E 3 ,„ 9, 2 5 .md'^.r, 7. , , 

7 = -[Y5(6^')^^^-1-01(^)'^^]- (9) 

Since the second r.-h.-s. term in Eq. (|9]) is negative, and its magnitude increases faster than 
the kinetic energy when the density increases, there exists a critical density Hc for mechanical 
collapse, {diildn)nc = 0, where the chemical potential /x is given by 

dE rl , 9s 2 2 ,m(i^,9 4, . . 

As shown in Fig. [H the chemical potential reaches the maximum at the critical density 
Uc = 2.04(/i^/m(i^)^, whereas the simple HF result implies ric = 9. 97 [fi^ / {md'^)]^ . Beyond 
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FIG. 1: Chemical potential /x as the function of the fermion density n, where the unit of fi is 
Ed = fi^ / {rm? d!^) and the unit of n is = \Ji?' / {md^y^ . The solid and dashed lines are the 
chemical potentials with and without the contribution from correlation energy respectively. 



this critical density ric, the compressibihty turns negative and the system is no longer stable 
against mechanical collapse. We should point out that the perturbation method is valid 
when the absolute value of the interaction energy is smaller than the kinetic energy, n < 
9 .58(71^ / md"^)^ as indicated by Eq. (Q. At = 2.04(/i^/m(i^)^, the perturbation still works, 
but higher order terms may be needed to produce more quantitatively accurate results. 
Nonetheless, the correlation energy effectively makes the system more vulnerable towards 
mechanical collapse. 

In the experiment on KRb molecules, the dipole moment for the singlet rovibrational 
ground state is d = 0.566(17) Debye [9[. Our calculation suggests that the critical density for 
mechanical collapse in this case is approximately ric — 9.23 x lO^^cm"^. In other references, 
this critical density is estimated to be 1.48 x lO^^cm^^ in the variational approach jlo| . 
4.89 X lO^^cm"^ from the zero-sound analysis jsj, and 4.80 x lO^^cm"^ in the collective- mode 
study Uj. The highest molecule density in the current experiment is about 0.3 x lO^^cm"^ 



12| when the temperature is at 1.4 times the Fermi temperature. It is hopeful that these 



theoretical results can be tested in future experiments with higher molecule density achieved. 
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III. TRAPPED DIPOLAR FERMI GAS 



The Hamiltonian of a dipolar Fermi gas in a harmonic trap is given by 



H ^ J dW(r)[__ + C/^jV;(r) 



+ ^ / I d'rd\'r{r)^p\r')V,,{r-r')^pir')^{r), 



(11) 



where the dipole-dipole interaction Vrfrf(r) is given by Vrfd(r) = {(P/r^)[l — 3(z^/r^)], and 
ip{r) is the fermion field operator. In the following we consider an axially symmetric trap 
with the trapping potential given by Uho — m{ujpx'^ + cOpi/^ + uj1z^)/2 and the trap aspect 
ratio given by A = uiz/^p- 

We consider the case that the scale of the trap is much larger than the interparticle 
distance, so the local density approximation (LDA) can be applied. The total energy E can 
be written as 

E = Ekin + Eho + Enar + 

where Ekim E^o, and Effar are kinetic, trap, and Hartree energies. The energy Es includes 
both the second-order Fock and correlation energies. Density functionals of these energies 
are given by 



Eho = / dhUho{v)n{v), 
EHar = \j Jd\d'r'n{r)V,,{r-r')n{r'), 
^2 j dhn^r) 



E, = -1.01- 



where n{r) is the fermion density with the number constraint N — J d^rn{r). At the ground 
state, the fermion density n{r) satisfies the energy-extreme condition, 

5N ^ 

from which the chemical potential n = dE/dN can be determined, 

^ = i^(67r2)lni(r) + ^(a;2 + y2 + AV) 
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- 2.36^n3(r). (13) 

In LDA, the total energy E satisfies a virial theorem which can be proved by the scahng 
transformation of the density [l^ 

n{r) r]^^°'n{r]r), (14) 

where rj and a are general numbers. Under this transformation, the total density and various 
parts of the total energy E transform as 



K 77 

kini 
hoi 
Har 1 



Es ^ r/^+i"E,. 

Close to ?7 = 1, to the first order 'm Srj = rj — 1, the changes in these quantities are given by 

5N = a6T]N, 
5 

SEkin = (2 + -a)5r]Ekin, 



5Eho = i~2 + a)6r]E, 



hoi 



SEnar = {3 + 2a)6r]EHar, 

SE, = {A+^-a)6r]E,. 
Therefore the change in the total energy is given by 

6E = 6Ekin + SEho + 6EHar + SEs, (15) 



and from the energy-extreme condition Eq. f[T2|) it also satisfies the equation 



= g^6N. (16) 



From the above equations, we have 



dE 5 
"^dN ^ i'2 + ^a)Ekin + i-2 + a)Eho 

+ i3 + 2a)EHar + i^+la)E„ (17) 



which leads to two independent equations 



2Ekin — 'iEho + ^Euar + 4-E's = 0, (18) 

5 7 BE 

-Ek,n + Eho + 2EHar + 'Es = N — . (19) 



Eqs. ( fT8l) and ( lT9ll together with identities 



dE 

d— = 2EHar + 4S„ 

od 
dE 

^pT. — 2Ehoi 

dujp 
dE 

fn— — = Eho — Ekin + Eg, 
am 

yield 

^dE _ 4^9E 

dd 3 dm ' 
1 

'^p^~ = E + i^m—, 
oojp 6 om 

dE A I dE 
N-— = -E + -m—, 
dN 3 9 dm' 

indicating the total energy E of the form 

E{m, Up, d, N) = hupNlgiN-^Cdd), (20) 

where Cdd is the dimensionless dipolar interaction strength, Cdd = md'^/aoh'^, ao = ^Jh/mcOp, 
and g{x) is a general functional of the variable x. 

In a trap, the system also suffers mechanical collapse if the total number of fermions 



exceeds a critical value. In the experiment on KRb polar molecules js], the optical trapping 
frequency for the polar molecule is about Up/ {2tt) = 109 Hz and the permanent electric dipole 
moments measured with stark spectroscopy is 0.566(17) Debye for the singlet rovibrational 
ground state, corresponding to dimensionless dipolar interaction strength Cdd = 0.2844. 
Our numerical results for the experimental condition show that the critical total number of 
polar molecules depends on the shape of the optical trap. As shown in Fig. when 
is fixed and A varies, the critical total number of polar molecules monotonically increases 
with the trap aspect ratio A, indicating that a more pronounced pancake shape of the trap 
can help stabilize the dipolar Fermi gas. Our results show that the critical molecule number 
for mechanical collapse is about 1.75 x 10^ for A = 1, slightly smaller than the HF result 
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FIG. 2: The critical total number A^^^ri of polar molecules for mechanical collapse versus the 
trap aspect ratio A for the singlet rovibrational ground state of KRb with fixed longitudinal trap 
frequency uJz/{2tt) = 109 Hz. 



2.15 X 10^ 



151 ]. implying the correlation effect further destabilizes a trapped dipolar Fermi 



gas. Again we would like to make a remark that our result about the correlation energy 
given by Eq. ^ is obtained in the second-order perturbation theory. Better results about 
the correlation energy are needed to estimate the critical molecule number for mechanical 
collapse with more accuracy. 



IV. CONCLUSION 



We studied the normal state of a dipolar Fermi gas beyond the HF approximation. For the 
homogeneous system, we obtained the ground-state energy including the correlation energy 
up to second order in the dipole interaction strength, which leads to a critical density for 
mechanical collapse smaller than that estimated in the HF approximation. For the trapped 
dipolar Fermi gas, we proposed a energy functional to describe this system based on LDA. 
We show that this energy functional satisfies a virial theorem. The instability of this system 
is also investigated. In both cases, the correlation energy makes the system more vulnerable 
toward mechanical collapse. 



9 



Acknowledgment 



We would like to thank T.-L. Ho, S. Yi, and H. Zhai for helpful discussions and other 
colleagues for constructive feedbacks. This work is supported by NSFC under Grant No. 
10974004, and by Chinese MOST under grant number 2006CB921402. 

REFERENCES 



[1] L. D. Carr, D. DeMille, R. V. Krems, and J. Ye, New. J. Phys. 11, 055049 (2009). 

[2] M. Baranov, Phys. Rep. 464, 71 (2008). 

[3] T. Miyakawa, T. Sogo, and H. Pu, Phys. Rev. A 77 061603 (2008). 

[4] L. You and M. Marinescu, Phys. Rev. A 60, 2324 (1999). 

[5] S. Ronen and J. L. Bohn, Phys. Rev. A 81, 033601 (2010). 

[6] W. Kohn and J. M. Luttinger, Phys. Rev. 118, 41 (1960). 

[7] T. D. Lee and C. N. Yang, Phys. Rev. 105, 1119 (1957). 

[8] W. H. Press, S. A. Tcukolsky, W. T. Vetterhng, and B. P. Flannery, Numerical Recipes in C 
(Cambridge, New York, 1986). 

[9] K. K. Ni, S. Ospelkaus, M. H. G. de Miranda et al, Science 322, 231(2008). 

[10] T. Sogo, L. He, T. Miyakawa, S. Yi, H. Lu, and H. Pu, New. J. Phys. 11, 055017 (2009). 

[11] C.-K. Chan, C. Wu, W.-C. Lcc, and S. Das Sarma, Phys. Rev. A 81, 023602 (2010). 

[12] K. K. Ni, S. Ospelkaus, D. Wang et al, Nature 464, 1324(2010). 

[13] K. Coral, B.-G. Englert, and K. Rzaiewski, Phys. Rev. A 63, 033606 (2001). 

[14] K. Coral, M. Brewczyk, and K. Rzaiewski, Phys. Rev. A 67, 025601 (2003). 

[15] J.-N. Zhang and S. Yi, Phys. Rev. A 80, 053614 (2009). 



10 



